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ABSTRACT: We discuss the possibility that heavy fermion superconductors involve 
odd-frequency triplet pairing. A key technical innovation in this discussion is a Majorana 
representation for the local moments which enables the nj = 1 constraint in the Kondo 
lattice to be handled without a Gutzwiller approximation. Our mean field theory for odd 
frequency pairing involves a condensation of local moments and conduction electrons, and 
is characterized by a spinor order parameter. This is a stable realization of odd frequency 
triplet pairing. It predicts a Fermi surface of gapless quasiparticles, whose spin and charge 
coherence factors vanish hnearly in energy. The unusual energy dependence of coherence 
factors leads to a T"^ NMR relaxation rate that coexists with a linear specific heat. A 
prediction of the theory is that a Korringa relaxation will fail to develop in heavy fermion 
superconductors, even in the limit of strong pair-breaking and severe gaplessness. 
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1. Introduction 

Heavy fermion superconductivity has attracted great interest in recent years as a candi- 
date for electronically mediated pairing. ^'"^ Six heavy fermion metals arc superconducting 
at room pressure: CeCu2Si2^, UBeis'^, UPt^^, URu2Si2^, UNi2Al^ and UPd2AlsJ 
These metals contain a dense array of magnetic rare earth or actinide ions that collec- 
tively participate in the formation of the superconducting state. Many properties set these 
systems apart from traditional superconductors. In particular, power laws in the specific 
heat, thermal conductivity, NMR relaxation rate and acoustic attenuation all point to the 
existence of gap nodes, and have been interpreted in terms of gap zeroes along lines of the 
Fermi surface. Furthermore, each of these superconductors appears to coexist with some 
measure of antiferromagnetic order. 

Existing phenomenological models of heavy fermion superconductivity treat it as a 
pairing process involving the pre-formed heavy f-quasipar tides. The strong repulsive 
interactions between these f-quasiparticles favor the development of nodes in the pair 
wavefunction, as suggested by the preponderance of power laws. Theoretical work on heavy 
fermion superconductivity has focussed largely on the possibility of momentum anisotropy 
in the gap function Aj^ as the origin of this node formation. The simplest candidates for 
gap functions with nodes are odd-parity triplet pairing, = — A_g, or even parity d-wave 
pairing. Two points appear to favor the latter possibility: 

o d-wave pairing is favored by the antiferromagnetic interactions that are characteristic 
to heavy fermion compounds. ^'■'^"^ 

o lines of gap zeroes inferred from many power-law properties of the condensed state, e.g 
T"^ dependence of the NMR relaxation rate, dependence of specific heat and ther- 
mal conductivity, effectively rule out triplet odd-parity pairing. Simple symmetry argu- 
ments show that odd-parity triplet would give rise to a gap vanishing at points, rather 
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than lines on the Fermi surface, in the presence of strong spin-orbit scattering. ■'^^'-'^^ 

Unfortunately, there are several observations that do not fit naturally into the d-wave 
scenario. One puzzling observation is the persistence of the NMR relaxation rate in 
heavy fermion superconductors with very large densities of gapless excitations and corre- 
spondingly "heavy" linear components to their specific heat.-*^^'-*^^ This gaplessness has been 
attributed to pair breaking by resonantly scattering off non-magnetic defects. Remark- 
ably, the constant density of quasiparticle states never appears as an observable Korringa 
NMR relaxation. In UPt^ and Ui^xThxBei^ for example, ■'^^'■'^^ the linear specific heat is 
of the same order as the normal phase value, yet there are more than two decades^^'^^ of 

spin relaxation. 

C 1 

7 -|- BT > linear term from pair breaking. 



T 



= ? -|- DT^ ^ no linear term from pair breaking? 



TiT 



(1.1) 



In conventional gapless superconductivity, spin coherence factors are unity at the 

Fermi energy. The robustness of the T NMR signal suggests vanishing spin coherence 
factors: a feature not easily accommodated by a conventional pairing hypothesis. 

The d-wave scenario is also unable to explain the isotropy of the H—T phase diagram of 
U Pt^}^'^^ U Pt^ has three separate low temperature flux phases that have been interpreted 
in terms of anisotropic pairing. There is a two-stage phase transition at zero field associated 
with the symmetry breaking effects of the weak heavy fermion antiferromagnetism, and 
perhaps also, a recently discovered incommensurate charge density wave.^^ The d-wave 
scenario supposes a gap function that transforms under a two-dimensional representation 
of the point-group: though this picture can account for the two stage transition, it predicts 
a two-phase flux lattice for all orientations of the applied field. 
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Finally of course, the d-wave pairing picture of heavy fermion superconductivity makes 
no reference to the close link between heavy fermion superconductivity and magnetism. 
Typically, the entropy associated with the superconducting phase 



is a significant proportion of the Rln2 entropy associated with the quenching of the low lying 
doublets: in this sense heavy fermion superconductivity is a spin ordering process, involving 
the magnetic, rather than the charge degrees of freedom of the f-electrons. Experimentally, 
each heavy fermion superconductor appears to have a coexistent antiferromagnetic order. 
In the recently discovered 1 — 2 — 3 compound UPd2Al^, an ordered moment of O.Sfig 
coexists^^ with the superconductivity. In URu2Si2, there is also evidence for a large 
moment-free order parameter that breaks time reversal and translation symmetries. ^^'^^ 
Unlike the well-known Chevrel phases, this moment shares the same magnetic degrees 
of freedom that are involved in pairing. It is rather difficult to account for coexistent 
magnetism and superconductivity in terms of two weakly coupled order parameters. 

These difficulties motivate us to reconsider the way in which heavy fermion supercon- 
ductors develop nodes in the pair wavefunction. Past analyses of heavy fermion supercon- 
ductivity have focussed on the spatial anisotropy. In this paper we explore a new avenue, 
examining the possibility of pair condensation into a state where the pair wavefunction 
has odd temporal parity. In this hypothetical state, pairing is retarded and the pair 

wavefunction contains a node in time.'^'^ 

Berezinskii^^ first pointed out that a general pairing hypothesis must consider the 
symmetry of the pair wavefunction under frequency inversion. Let us denote the pair 
wavefunction 



Here (. . .) denotes the time-ordered expectation value, and we use a four-vector notation 
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(1.2) 



K = (k, w). Since the Fermi operators anticommute, the pair wavefunction satisfies 

F{k) = (1.3) 

where [F^]q,^ = [F]f3a denotes the wavefunction with spin indices transposed. Now if we 
assume that the state breaks neither time reversal symmetry, nor spatial parity, then the 
pair wavefunction must have distinct spatial, temporal and spin parity. Let (S, P, T = 
±1) be the parities of the pair wavefunction under the interchange of spin, space or time 
coordinates respectively, i.e. 

'SF'^(k,a;) 

F(k,a;) = I PF(-k,a;) (S, P, T = ±1) (1.4) 

. TF(k, -u) 

then the total antisymmetry of the pair wavefunction implies that the combined product 
of all three parities must equal —1: 

SPT = -1 (1.5) 

Superconductors with an antisymmetric spin wavefunction, S = —1 are "singlet" super- 
conductors, 

F{k) = ia2Fs{K) {F = -FT) (1.6) 

whereas superconductors with a symmetric spin wavefunction, S — +1 are "triplet" 

F{k) = ia2l- Ftin) {F = +F^) (1.7) 

In conventional superconductivity, T = +1, so that the spatial parity of singlet and triplet 

states is even and odd respectively. Berezinskii has argued that symmetry also permits 

the possibility of odd frequency pairing where 

+ 1 (triplet {S,P,T}^ {+,+,-}) 

T=-l P={ (1.8) 
-1 (singlet {5, P,T} = {-,-,-}) 
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Odd frequency, even parity triplet pairing was first considered by Berezinskii in the context 
of He — 3. A renaissance of interest in these types of states has been prompted by the work 
of Balatsky and Abrahams, who are the first to discuss the possibihty of odd frequency, 
odd parity singlet pairing. 

Historically, odd frequency pairing has not enjoyed a great deal of attention. One 
reason for this lack of attention is that the simplest odd frequency paired state is luistaljlc , 
with a negative Meissner phase stiffness. For example, in the s-wave triplet state, (S, 
P= +1) the momentum dependence of the gap function vanishes, and the London Kernel 
is formally identical to s-wave singlet pairing 



Since an odd gap function where An = —A-n must also satisfy the analyticity requirement 
A^j = Al^, this implies A^j is purely imaginary. Thus A^ < and the stiffness is negative. 
Loosely interpreted, this negative stiffness suggests the microscopic phase of the order 
parameter likes to "coil up", breaking translation symmetry and assuming a staggered 
configuration. 

We shall argue that the Kondo effect between a conduction sea and local moments 
in heavy fermion metals provides an ideal source of retarded scattering for odd frequency 
pairing. In the normal state, this retardation generates resonant bound states between the 
conduction electrons and local moments, quenching the moments and forming the heavy 
quasiparticles. In the superconducting state, the resonant Kondo scattering acquires a 
pairing component that results in even parity, odd frequency triplet pairing of the conduc- 
tion electrons. This state develops a phase stiffness by the simultaneous condensation of 
the local moments and the conduction electron pair degrees of freedom. The equal time 
order parameter is a matrix correlating these two degrees of freedom 




n 



'n 




(1.9) 




{a, P= 1,2,3) 



(1.10) 
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Here S[xj] denotes the local moment spin at site j; t[x] is the conduction electron "isospin", 
whose z component describes the number density, and transverse components describe the 
pairing 

T3 =i^{p{x) - 1) 

t t 

r+{x) =^/;^(x)V'| {x) 

The quantity g defines the magnitude of the order parameter. is an orthogonal matrix 
whose rows define an orthogonal triad of unit vectors dx, 

( di{x) \ 

M{^) = d2{x) . (1.12) 
V dsix) I 

Stabilization of the odd-frequency paired state is achieved by staggering this order pa- 
rameter: in a simple model, di and d2 are staggered commensurately with the lattice. A 
"composite" order parameter of this form has been recently suggested, in the context of 
the two-channel Kondo model, as an order parameter for odd frequency pairing, by Emery 
and Kivelson."^^ 

Within our theory, the microscopic manifestation of this type of pairing is an anomalous 
self-energy in the triplet channel with a pole at zero frequency 

A{k) = ia2(k (4 = [dl + id2] ■ a) (1.13) 

A spinless component of the conduction electron band decouples from this singular pairing 
field, giving rise to surfaces of gapless excitations. Spin and charge coherence factors of 
these quasiparticles vanish linearly with the energy on the Fermi surface, 

(k'|Pk_k'|k> 



oca; {u; = E^, + E^) (1.14) 



creating the unusual circumstance where a fiow of quasiparticles transmits heat without 
passage of charge or spin. These unusual coherence factors lead to power laws in the 
nuclear magnetic relaxation 

^ocT3 (1.15) 



that coexist with a hnear specific heat capacity. 

An essential part of our analysis is a second quantized description of the local moments 
that avoids constraints. Existing treatments factorize the spin variable in terms of spin 
1/2 fermions. 

S = f'^a[^Uf0 (1-16) 

This approach requires a constraint n^- = 1 to impose the condition S = -^j which is 
the origin of additional complications. In many practical applications, the constraint is 
weakened, imposing it at the mean field or the Gaussian level of approximation. Here, we 
employ real, or "Majorana" fermions to represent spins. Perhaps the most famous example 
of a Majorana fermion is a single Pauli spin operator. Recall that 

so the fermions 

-L.- (1.17) 

also satisfy a canonical anticommutation algebra = 5^5. Indeed, for any such 

triplet of Majorana fermions, it follows that the "Spin operators" 

S=~ffxff (1.18) 

simultaneously satisfy both the spin algebra, and the constraint S"^ = 3/4. The general- 
ization of this result to a lattice of spins, employs a vector of Majorana fermions ry^ = 
(a = 1, 2, 3) for each site j where, 

{vt,V^j} = SijS-^ (1.19) 
from which the spin operator at each site is constructed 

= ~Vj X ffj (1.20) 
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These spin operators behave as independent spin 1/2 operators. 

Our Majorana representation of spins provides a natural lattice generalization of an- 
ticommuting Pauli operators. On a lattice, we may represent Majorana fermions in terms 
of a set of N/2 independent complex fermions that span half the Brillouin zone: 

fij = -^ E {^k^^'^-^^+^V^'^-^^} (1-21) 

kGl/2B.Z. 

These complex fermions obey canonical commutation relations {rj^^i]^^} = 6"'^6,-',-',; op- 

k k' 

posite halves of the Brillouin zone are related as complex conjugates: = v"" Since 

the Fock space is spanned by 3N/2 complex Fermi operators, it is 2^/^ times larger than 
a Hilbert space of N commuting spin 1/2 operators. The spin algebra and the condi- 
tion S — 1/2 are satisfied between all states of the Fock space, thus the anticommuting 
representation replicates the spin Hilbert space times. "^^ We may then represent the 
partition function of an electronic system containing N spins as an unconstrained trace 
over the independent Fermi fields 

^ = K^^'^'^^n ^ -'-rf, X rj,) (1.22) 

where the formal normalization factor associated with the replication of states has been 
added. 

Our basic model for a heavy fermion system is an 5 = 1/2 Kondo lattice with a single 
band interacting with one local f-moment Sj in each unit cell. In a real Kondo lattice, 
the local moments are strongly spin-orbit coupled into a state of definite J. We shall 
assume that the low lying spin excitations are described by a Kramers doublet, where a 
low energy S = 1/2 Kondo model becomes more appropriate. For simplicity, we shall 
ignore the anisotropics that are necessarily present in a real heavy fermion system. In our 
model, the latent superconducting pairing is driven by the on site Kondo interactions, and 
the state that forms exhibits a coexistence of magnetism and superconductivity. 
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The simplified isotropic Kondo lattice model that we shall use is then written 

H = Hc + J2Hint[j] (1-23) 
3 

where 

^^ = E^kV'^kV'k (1-24) 
k 

describes the conduction band, and tp^^ = (V'^^-|-) V'^^p is a conduction electron spinor. 
The exchange interaction at each site j is written in a tight binding representation as 

HintiJ] = J {'iP^ja^ap^jf3) ' (1.25) 
When written in terms of the Majorana fermions, this term becomes 

HinM = -^^"^A^ ■ V'i ^ (^^:| ) (1.26) 

where we have used the result ia-lff xif) = [^7 • i?]^ ^ | to simplify the interaction, absorbing 
the bilinear term as a redefinition of the chemical potential. This simple form of the 
interaction can be rewritten in a suggestive form, by defining the composite spinor operator 

The Kondo interaction is then the "square" of this operator: 

HinAJ] = -J^h^j (1-28) 

suggesting that in the lattice, we should consider the possibility of states where the local 
moment and electron spins condense together to develop a vacuum expectation value of 
this spinor quantity 



This order parameter transforms as a spin 1/2 object, so changes in its sign correspond to 
physical rotations of the condensate by 27r. Defects in the spinor field are then disclination 
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lines or "Z2 vortices" , around which the phase of the spinor order parameter changes by 
TT (Fig. 1). The gauge equivalent integral of the vector potential around a Z2 defect is 

^jA-dx = 7r (1.30) 

so the flux quantum of a "charge e" spinor is the same as a charge "2e" scalar: 

^0= I A-dx=-TT=^ (1.31) 
J e 2e 

In our model there is a microscopic "Z2" gauge symmetry 

fjj ^ ±r/j (1.32) 

Vj transforms in the same way as the Majorana fermions and thus its phase is defined 
to within ±7r. Physical quantities involve the combinations of the square of Vj at each 
site and are Z2 invariants. The three Z2 invariant quantities that this defines are just the 
components of the matrix M.{x) 

{p{x)S{x)) =gV{x)'^aV{x) 
{'4!\{x)^\s{x)) =gV^{x)ia2BV{x) (1.33) 

Under a phase change of tt in the spinor field, the axes of the composite order parameter 
rotate through 27r. 

The outline of this paper is as follows: 

2. Development of a path integral formulation of the Kondo lattice, demonstrating how 
the simplest decoupling procedure leads to an odd frequency paired state. 

3. Discussion of the quasiparticle excitations and coherence factors. 

4. Calculation of the mean field thermodynamics in this paired state. 

5. Computation of the Meissner stiffness of this phase, and form of the Landau-Ginzburg 
theory. 
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6. Effect of vanishing coherence factors on local magnetic and charge response. 

7. Interplay with magnetism. 

8. Critique and discussion: possible application to the theory of heavy fermion supercon- 
ductivity. 

Certain formal arguments, not pertinent to the main flow of ideas have been reproduced in 
the the appendices. In appendix A, we show how the Majorana representation is related 
to the Abrikosov pseudofermion representation. In appendix B we give some examples of 
the application of the Majorana representation to simple spin models, showing the relation 
to the Jordan Wigner transformation in the one dimensional Heisenberg model. 

2. Path Integral Representation of the Kondo Lattice Model 

To develop a "toy model" for the odd-paired state, we focus our attention on a stripped- 
down Kondo lattice model, with the Hamiltonian described in (1.23) to (1-26) 

^ = E ^k^^k^k - i E V'^- ■ (2-1) 

k j 

We have suppressed both the momentum dependence and anisotropy of the coupling. In a 
real heavy fermion system, we envisage that the spin indices would refer to the conserved 
pseudospin indices of the low lying Kramers doublets. 

To illustrate the calculations in this section, we shall use Feynman diagrams, as shown 
in Fig. 2. The bare propagator for the conduction electrons is represented by a solid arrow, 
the bare propagator for the Majorana fermions by a dashed line, without an arrow. 

— — — — =('?''(«)'?'''(-«))<. = *''''(—') 

The product form of the exchange interaction (1.26) clearly suggests a decoupling in terms 
of the spinor variable 

^i = (l/:i)=— ([^•^iiv'i) (2.3) 
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corresponding to the bound state of an electron and a local moment. 



With this point in mind, we now write the partition function as a path integral, Z = 
Jj>e~ -^0 ■^W'^'^ where 

AT) = EV^k^^-V^k+ E V^'^drV^ + Hc + J^Hintij] (2.4) 
k ke^BZ 3 

and we have factorized the interaction in terms of a fluctuating two-component spinor 

v\j = {v*,vp 

Hintb] = ^"^ji^ ■ Vj)Vj + V'fjia ■ fjj)^j + 2\Vj\^/J (2.5) 

For later purposes, it is particularly useful for us to introduce a Balian Werthammer four- 
spinor notation, defining 



/ ^3} \ 



(2.6) 



The lower two entries of each spinor are the time reversed pairs of the upper two entries. 
In terms of these spinors, the conduction electron Hamiltonian is 



keiB.z. 



where 



(ti, T2, T3) = 



1 



—I 



-1 



(2.7) 



denotes the triplet of isospin operators. Note the definition of T2. The factorized interaction 
can also be written 



Hint\j] = 2 [*^(^ • + Vtj (a • fjj)^j\ + Vjhj/J 
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(2.8) 



where, 



a 



g_ = a ® \ — 



(2.9) 



denotes the spin operator in the Bahan Werthammer notation. 

We are particularly interested in expanding around static mean field configurations 
where the amplitude of Vn is constant, 



This choice of mean field theory is equivalent to a resummation of the interaction lines in the 
pairing channel between the conduction and Majorana fermions, leading to a saddle point 
condition for the anomalous average of Majorana and conduction electrons, as illustrated 
diagrammatically in Fig. 3(i). The development of this anomalous average leads to self- 
energy insertions in the conduction electron lines (Fig. 3(ii)). From this diagram, it is 
evident that the conduction electron self energies are bilinear forms in the spinor Vj, and 
are hence mvariant under the Zi gauge symmetry. 

We can actually find a class of degenerate mean field solutions by arbitrarily reversing 
the sign of Vj — > TnjVj (rrij = ±) at any site. For each choice of sign, there are 

2N/2 

equivalent ways of choosing the independent Majorana creation operators in momentum 
space, {v^ _^ = ■^j^), there are thus 2-^"-^/^ = 2-^/^ independent degenerate saddle point 
solutions for each static solution {Vj}. Each saddle point is physically identical, so we may 
absorb the 

2N/2 

normalization in the partition function by restricting our attention to one 
representative saddle point 



In this way, we fix the gauge for the local Z2 invariance. 

On a bipartite lattice, it is convenient to carry out a gauge transformation that moves 






(2.10) 




(2.11) 
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the origin of momentum space to the Brillouin zone center, by defining 



9j=Q-Rj (2-12) 



11 ^ IT, ~ 



The conduction electron Hamiltonian can then be written 



Hc= *"^k[%-^k^3]*k (2.13) 



keiB.z. 



where 



k~2^k-Q/2 "k+Q/2J 

_1 , 



(2.14) 



and we have suppressed the tildes on the electron operators. For a simple bipartite tight- 
binding lattice, taking Q = (7r,7r,7r), then = n and = ^^_q^2' ^^^^ ^^^^ choice of 
gauge, the conduction electron kinetic energy is manifestly particle-hole symmetric. 

The lowest energy mean field solution is obtained for a spatially uniform Vj = Vo- With 
this choice, the admixture between conduction electrons and local moments is described 
by the mean field Hamiltonian 

i^mzx^ E [*V^-^k)Va + vt,(r7tj^.a)vE'^] (2.15) 
ke|BZ 

To gain further insight into the nature of this state, let us "integrate out" the Majorana 
fermion degrees of freedom. Within a diagrammatic approach, this corresponds to intro- 
ducing a self-energy into the conduction electron propagators, as shown in Fig. 3(ii). The 
effective action for the conduction electrons contains this self energy: 

^eff= E («) [^- eg -//£T3-E («)]*(«) (2.16) 
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S(k) describes the resonant scattering through zero energy spin states, and is given by 



(2.17) 



where the matrix is formed from the square of the spinor order parameter 
The general expansion of A4 is 

M = \v^\l + daba^®A 

where ^a^fj] = [j](T'^ <^t^ Z[j]. The columns of the matrix d define a triad of orthogonal 
unit vectors [d^]a = d^x^ — 1? 2, 3) 



dj + idj =Zj ia2crzj 



(2.18) 



that set the orientation of the order parameter in spin space. The resulting conduction 
electron self-energy is proportional to 1 /a; 



F2 

E(«) = —P 



(2.19) 



where the projection operator 



(E = V) 



(2.20) 



The anisotropic component of the self-energy S'^'^(k;) = ^['P — |] contains "anomalous" 



components, and may be written 



{B,{u) ■ a)n + (At(^) . a)^ + (A(c^) • af- 



(2.21) 



We interpret the quantities 

B{uj) =A{u)d^ 

A{u}) =A{uj){d^ + id^) 



(2.22) 
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as resonant exchange and triplet pairing fields, respectively. This is a realization of odd 
frequency pairing in the triplet channel. Unlike earlier realizations of odd frequency triplet 
pairing, ^9,30 ^j^g g^^p function diverges at zero frequency. Such resonant contributions to 
the self-energy are well known within mean field treatments of the Kondo lattice, but here 
the resonant scattering acquires additional, anisotropic pairing terms associated with the 
pair condensate. 

The projective form of the pairing self-energy means that not all components of the 
conduction fluid experience the resonant scattering. This is more easily seen by decom- 
posing the conduction electron operators into four Majorana components 

i^^^^[^l + i^^-a]zo (2.23) 

The "vector" components of the conduction electron are projected out by the operator V. 
By substituting into the mean field Hamiltonian, we see that only these components of the 
conduction electrons couple to the resonant scattering potential 

Hrmx= E -iV[P^-rf^-R.C.] (2.24) 
kG^B.Z. 

The zeroth component ^'o = (1 — £)^ does not couple directly to the resonant scattering 
potential, leading to a gapless quasiparticle mode. 

We may gain insight into the meaning of the order parameter matrix Al by making 
the identification 

Vj = -'^^-ffji^j (2.25) 

or 

Vj = -^a-ffj^j (2.26) 

inside the path integral. The expectation value of the matrix A4 is then the irreducible 
part of the corresponding product of operators. In particular, 

(Vj^r-a'Vj) = ^(*/r«(a • r,j)A^ ■ Vj)^j)l (2-27) 

17 



where "I" denotes the irreducible part. Using the identity {a-r]j)a{a-r]j) = [—if)jXf]j—a/2] 



it follows that 



Vj) = j''{T-[xj]S%]) 



(2.28) 



where S{xj) = Sj is the local moment at site j and 



(2.29) 



is the conduction electron "isospin". This order parameter represents a bound-state be- 
tween the local moments and the conduction electron charge and pair degrees of freedom. 
The composite order parameter 



represents the development of a joint correlation between the conduction electron singlet 
pair density and the local moment spin density. Clearly, this state breaks (i) electron 
gauge symmetry, (ii) spin rotation symmetry and (iii) time reversal symmetry. Despite 
these features it does not necessarily follow that the state formed has either an ordered 
moment, or an equal time pairing field. 

3. Excitation spectrum and quasiparticles 

Let us now examine the nature of the excitation spectrum in this odd paired state. Let 
us begin by rewriting the mean field Hamiltonian in terms of the Majorana components 
(2.23), then 



dc{x) = (P{x) — i(P{x) 



{^P^{x)iP^{x)S{x)) 



(2.30) 




(3.1) 



V&\B.Z. 



where 



(3.2) 
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is the total charge operator. Let us consider the special case of = 0, when the spectrum 
of the zeroth component remains unrenormalized. In this special case, the Hamiltonian 
can be written in terms of quasiparticle operators as follows 

H =Ho + Hg 

keiB.Z. (3.3) 

^9= E %«'^ka«ka 
k, a=(l,2,3) 

The first term describes a gapless "Majorana" conduction band that spans the half Bril- 
louin zone where > 0. The second term describes a gapped band with excitation energies 



This band spans the entire Brillouin zone, since it incorporates three Majorana conduction 
and three Majorana spin fermions. The basic character of the quasiparticle spectrum is 
unchanged when we consider finite deviations from particle- hole symmetry //^ 7^ 0. There 
are two important features (Fig. 4): 

— * — * 

1. A three-fold degenerate gapful excitation centered around k = —Q/2. In the vicinity 
of the gap 



^9 "I ~f 



k Q/2, a=(l,2,3) 

where ~ V'^ /D and m* = -^rn and m is the conduction band mass at the band 
edge. Just above the gap, the quasiparticles have almost no conduction character: 
correlation functions of the local moments are then determined through the relation 

kel/2B.Z. 

— * 

2. A neutral Majorana band, located around k = 0. 
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To gain more insight into these excitations, let us consider the conduction electron 
propagator 

g-\oj) = [a; - - - (3.7) 
Choosing Zq = ^ Q ^ ' ^^en [d^, d^, d"^) = {x, y, z), so that cZ^^ = ^^j^ and 

S(«) = — P; P = 7[3(l)-a"®T'^] (3.8) 
It is useful to define "up" and "down" spin projection operators 



=^[1 + ^3®T3] 



(3.9) 



2' 

We can use these operators to project out the "up" and "down" electron propagators. 
{GaPa = gPa)) 

G^=[{u-~e^-A,j) + HTs + AujTi]-^ 

A^ = A(a;) (3.10) 

It is also useful to evaluate the determinants 

det[GT\k,u)] =[{uj - ir - A^f - - AJ] 

(3.11) 

det[G-i(k,c^)] =[{uj - - 2A^f - 

Zeroes of these functions determine the quasiparticle excitation energies uj^^: 
det[G'-l(k,a;j^J] = 0(Fig. 5). 

The "down" propagator contains no pairing terms, and describes a gapful band of 
quasiparticles with excitation energies 



= ± y ' + y (3-12) 

This spectrum closely resembles the large N solution to the particle-hole symmetric Kondo 
model, with a hybridization gap 2A^.^^~'^^ 
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The "up" electron propagator describes a band of odd-frequency paired electrons. The 
poles of this propagator 



(3.13) 



at a; = a;^ are determined by the cubic equation det[G| ""^(k, a;j^)] = (see (3.11)). Solving 
for the conduction electron energy as a function of u>^, e^^ = e±(a;^), we find 



(3.14) 



which defines two branches of the "up" qiiasiparticle excitation spectrum. The (-) branch 
is gapful with a gap 2Ag ~ ~ ' ^^^^^ ^ conduction electron half- 

bandwidth. The (-I-) branch is gapless, corresponding to the Majorana component of 
the conduction sea that decouples from the resonant scattering center. The quasiparticle 
density of states corresponding to these two branches is 



de±{u!) 



1 + 



'-^co 



It 



At low energies 



giving an enhanced density of states A^^(O) = 2^ where Zq 
surface. 



1 + K 



(3.15) 



(3.16) 
at the Fermi 



Let us explicitly construct these gapless quasiparticles. It is convenient to split the 
Hamiltonian into "up" and "down" spin parts H = + H^, where 



k 
1 



(3.17) 



^ki=7i[^£ + ^^l 
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describes the hybridized band of unpaired "down" electrons, and 



V_ 
^/2 



V_ 
^/2 



_V_ 

_v_ 





(3.18) 



describes the paired "up" electrons. In terms of the quasiparticle operators 



ka ka ka 



(3.19) 



where only positive energies enter into the Hamiltonian. Here a = 0, 3 denotes the gapless 
and gapful excitation branch respectively. The quasiparticle operators for the gapless "up" 
electrons can be constructed from a generalized Bogoliubov transformation 



The eigenvector 0r containing the Bogoliubov coefficients satisfies 



(3.20) 



Mk = ^ko'^k 



(3.21) 



Eliminating Zr and substituting back into (3.21) then gives 



^f'fta;^^)(;k)=o 



(3.22) 



where ^{k) is taken from (3.10). Diagonalizing this eigenvalue equation gives = 
^^('^ko^' ^k = "^^(^ko)' ^^^^^ 



1 + 



(3.23) 
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where the energies are given by (3.14). = Z-\-{u^^) takes the form given in (3.15). 
Energy, rather than momentum-dependent Bogohubov coefficients are a characteristic of 
odd-frequency pairing. 

Let us now consider the charge and spin coherence factors of these gapless excitations. 
Within the gapless "up" band only the conduction charge and spin operators contain 
diagonal matrix elements. Suppose we attempt to excite quasiparticles out of the ground- 
state, by coupling to a charge, or spin excitation; the relevant diagonal matrix elements 
are 

(k-| { ^1 } |k+) = (k-|{V"^i^-^V'k+^ - V'_k-tV'"^_k+^}|k+) (3.24) 

where = k ± g/2 and |k^) = ot-j_|0) denotes the state with one quasiparticle added to 
the mean field ground-state. In terms of the Bogoliubov coefficients this is 

(k-, I I I |k+) = ^Z^^Z^_ [u+u- - V+V-] (3.25) 

where u± = u+{uj^^) and v± = v+{u!^^) are the Bogoliubov factors in the gapless band 
(-I-). On the Fermi surface, since u — v = this coherence factor vanishes . Away from 
the Fermi surface, the spin/charge coherence factor grows linearly with energy 

(^~\{al]\^^) = H^+^^-){v^) (3-26) 

(In the special case of particle-hole symmetry, these coherence factors vanish throughout 
the gap.) In a similar fashion, we may examine quasiparticle components of the charge 
and spin, given by 



2 2 



'^jt:^ K « 
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From these results, we conclude that there is no way to couple via charge or spin probes 
to the quasiparticle excitations at the Fermi surface. These gapless excitations are devoid 
of charge, or spin quantum numbers on the Fermi surface. This dramatic effect is a direct 
consequence of the resonant pairing and the pole in the gap function. So long as this 
pole is maintained, the coherence factors will identically vanish on the Fermi surface. Note 
however, that these quasiparticles can still carry entropy, and in this sense can be regarded 
as thermal quasiparticles. 

It is particularly instructive to examine the local conduction electron propagator and 
the pair wavefunction in this simple mean field theory. The local propagator for the paired 
"up" electrons is 

rD 

G|(a;) = ^G|(k,a;) = p / deG^{e,uj) (3.28) 
— * «/ — D 



Carrying out the integral over the conduction electron energies we find 

1 {H>Ag) 



— Im 
7rp 




1 + sgn(a;) 



A(w)ti + ^T3 



(3.29) 



Loosely speaking, the electrons are normal outside the gap region and become paired at 
energies less than the gap A^. The spectral function can be rewritten in terms of the 
energy dependent Bogoliubov coefficients derived in (3.23) 



-Im 



7rp 



G-^{u> — iS) 



2 



(3.30) 



where the coefficients u^j and v^i are evaluated in the gapless band. 



Finally, we may construct the pair wavefunction from the off-diagonal components of 
this spectral function 

A(a;)sgn(a;) 

r 

(3.31) 



{ip^{x,iu;n)ip^ix,-iuJn)) =^ [ —-. — ^— 

2 J_J) 7T ILOn - 



= - liOnP 



^ VA(a;)^ + /x^ 
^ duj 1 A(a;) 



TT [a;2 + a;2] yA(ZJpT7? 
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thereby explicitly displaying the odd-frequency character of the pairing. 



4. Mean Field Thermodynamics 

Next, we discuss the mean field thermodynamics. The mean field free energy per site 
is written in terms of the conduction electron propagators as 



= ^-|j]in{det[^-l(«)]} 



(4.1) 



where the determinant can be expanded in terms of the "up" and "down" components of 
the propagator (3.10) 



det 



g-^in)] = det[G-^{K)]det[G-^{K)] 



(4.2) 



Differentiating with respect to the order parameter, yields the mean field equation 

2{iujn - eg - 2A„) 



J 2 iu!r,\i 



+ 



■}- 







(4.3) 



where the denominators in these equations are given in (3.11). For a constant conduction 
electron density of states p, we may replace 



(4.4) 



where the contour integral proceeds clockwise around the branch cut in the function 

' z - D ' 



-D-z 



(4.5) 



The energy integrals can then be performed by closing the contour around the poles in the 
Green functions, which are located at 



e|a(t^) =u; - + a sgn(a))^A2 + //2 



(« = ±) 



(4.6) 
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for the up and down electrons respectively. Carrying out the complex integral then gives 

where F(joi{z) = pQ[eaaiz)]/7r and we have used the notation = A{iujn)- Carrying out 
the Matsubara sums then yields 

where a;"*" = uj + i6. The functions F^^(u'^) = ImFa-a{uj~^) count the number of up and 
down excitation branches at frequency u). Ignoring the small differences between the up 

and down spin excitation gaps, 

r eiD-\aj\) (a=T« = +) 

FU^+) = p{ (4.9) 
[ 0{D - \uj\) - e{Ag - \uj\) (otherwise) 

for the gapless and gapful branches, respectively. This simple mean field theory then gives 
rise to a phase transition at a temperature 



Dexp 



= Tx exp 



(3/2 Jp) J - (6Jp). 
without the formation of an intermediate heavy fermion phase, where Tj^ — Dexp[— -^j^] 
is the single ion Kondo temperature. At = 0, the gapless excitation branch of the 
spectrum does not contribute to mean field equation. At finite the gapless branch 
develops a small linear coherence factor, and we see that this has the effect of suppressing 
the transition temperature. For all values of /j, however, the form of the mean field Ag(T) 
quite closely resembles that of a singlet BCS superconductor. 

The precise relation between the single ion Kondo temperature and Tc is not reliably 
predicted by the mean field theory. Our path integral approach amounts to a "Hartree" 
decoupling of the interaction. Had we chosen a more conventional diagrammatic approach, 
carrying out a "Hartree-Fock" decoupling of the original Hamiltonian, writing 

(#i)=-^^ ({/:{) (4-10) 
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so that 




(4.11) 



then the mean field Hamiltonian would have become 



Hintij] =^ i^^ji^ ■ Vj)Vj + V^jia ■ ffj)ijj + 3\V//2J 



(4.12) 



so that for this scheme, T^^^ = T^, to logarithmic accuracy. The path integral approach 
recovers the "Fock" contributions to the pairing as a leading order component of the RPA 
fluctuation corrections to the mean field transition temperature (Fig. 6). Since fluctuation 
effects will suppress Tc in either scheme, the particular choice of mean field theory is 
somewhat arbitrary, and will not matter at the next level of approximation. 

If we take our mean field theory literally, we see that from the point of view of the orig- 
inal conduction band, the transition into the odd frequency state can occur for arbitrarily 
weak coupling constant, taking advantage of the Kondo effect to produce a logarithmic 
divergence in the pairing channel. Significant pair breaking effects will of course come from 
the fluctuations. In the one impurity problem, there are infrared divergences in the Gaus- 
sian fluctuations that suppress the mean field transition temperature to zero. However, 
(unlike the large N approach) in the lattice model there is no continuous gauge symme- 
try so the development of a gap in the spectrum will cut off the one-impurity infrared 
divergences, preserving a finite temperature transition. 

Finally, we note that expanding the Free energy at low temperatures yields a linear 
specific heat proportional to the density of states in the gapless band 



Depending on the ratio (fx/D), this linear specific heat can range from a value characteristic 
of the free conduction band, {jJi/D) ~ to a value more characteristic of a heavy fermion 




(4.13) 



metal ~ 



^ for {^i/D) - 1 (Fig. 7). 
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5. Rigidity of the odd frequency paired state 

One of the key issues associated with odd frequency pairing, is whether it leads to a 
real superconducting Meissner eflFect. Past attempts to construct an odd frequency paired 
state within an Eliashberg formalism, have experienced difficulty in producing a state with 
a positive phase stiffness"^^ and a finite London penetration depth. To begin our discussion, 
we first discuss the form of the long-wavelength effective action. 

5. (a) Long-wavelength action 

In general, the mean field free energy will depend on gradients of the order parameter 
field, and the form of the applied vector potential. Let us now consider slow deformations 
of the order parameter 

z{x) = g{x) (^l^ (5.1) 

where 

— * 

is an SU (2) rotation matrix. The rate of rotation is given by V^r = giv, where 

uj = g~^Vg = iuj^s^ (5-3) 

is decomposed in terms of its components (jJx{x) along the principle axes dx of the order 
parameter. The leading quadratic terms in the gradient expansion of the Free energy about 
the uniform mean field theory are then 

A=l,3 

Here the stiffnesses px for slow twists about each principle axis are analogous to the 
moments of inertia of a top. 

To include the effects of an external magnetic field, we introduce a finite vector po- 
tential by an appeal to gauge invariance. Our original model is gauge invariant under the 
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transformation 

z{x) ^e*'^(^)^(a;) (5.5) 
Aix) ^l(a;)+ (^)v0(a;) 

so that 

g{x) ^^(a;)e*'^(^)^3 (5.6) 

This means that the long wavelength action must be invariant under the transformation 

LO^ — >Ld^ + 2V(/)(x) 

(5.7) 



-A(x) — '^Mx) + S/(f){x) 



In other words, a uniform vector potential A is equivalent to a uniform rotation rate —^A 
about the axis. The gauge invariant form of the Free energy is then 



F=^ I d^x 



Pm<^l + Ps{i:^3-^^f + —] (5.8) 
ti Ho J 

In terms of the vector h = d^ = z^az, this action can also be written 

F=- d^x pmi^hf + p,((^3 - -rAf + — (5.9) 

Z J I tl /Iq -I 

The term ps is the Meissner stiflFness of the superconductor, whereas the term pm can be 
regarded as a "spin stiffness" of the triplet paired state. It is at first surprising that a 
charge le spinor order parameter can give rise to a charge 2e coupling between gradients 
of the phase and the vector potential. We can resolve this apparent paradox by noting 
that that z{x) is also a spin 1/2 object, thus a phase change 

OJ ^ \0 

corresponds to a rotation through 9 = 2(p. A rotation through 27r leads to a sign change in 
z{x), and in physical configurations, z{x) must be continuous up to a sign ±1. The gauge 

— * 

invariant coupling between and A is then 
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and hence the couphng between physical rotations and the vector potential is a charge 
2e coupling, as in conventional superconductivity. Note finally, that if we include spin 
anisotropy into the original Hamiltonian, then this will tend to align the order parameter, 
for example, through the inclusion of a term of the form 



^ Pm f 



(fix{ds ■ zf (5.10) 

On length scales I > Iq, the system behaves as a conventional Landau- Ginzburg theory. 

5. (b) Computation of the Meissner Stiffness 

To compute the Meissner and spin stiffnesses in this gradient expansion, we consider 
a configuration with a uniform rotation about the principle axes da 



[S- ^ ^) (5.11) 



We may absorb this uniform rotation into a gauge transformation of the the conduction 
electrons through the replacement of the conduction electron kinetic energy by 

1 (5-12) 



Here we have expanded the kinetic energy to quadratic order in the twist. The effect of 
the twist in the phase can then be included into the electronic Green function by 

^ -g~\K) + h{K) (5.13) 

The Free energy of the system in the presence of a uniform twist can be calculated from 
the trace of the conduction electron propagator, as follows 

T 



F[coa] - F[0] = - -J2Tr[ln[-g-\K) + h{K)] - ln[-g-\K)]] 
= --J]Tv{/n[l-^(Ac)MAc)]} 
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Expanding the logarithm to quadratic order gives 



F[uj] - F[0] 



^"=«E(V)^ ^^^^ 



k 



(5.15) 



The first term in p'^ can be integrated by parts to yield 



(5.16) 



for the stiffness about the da axis. For our simple model, p = p = pm, and in the special 
case where // = all three stiffnesses are equal. We shall explicitly focus on the stiffness 
about the h = axis, which is associated with the London Kernel 

„2 



(5.17) 



Q^jj, =Ae^p2,5fj,y 

We can separate the trace into "up" and " down" components. The "down" component is 
unpaired and explicitly vanishes. The "up" component gives 



P3 



=4e 



2A: 



(5.18) 



£[(ic.„-A„-6^^)2-A2-^2]2 

In a conventional BCS theory, it is sufficient to impose a low energy cutoff on the 
frequency sum, 

\i^n\ < A 

after which, the conduction electron band-width can be taken to infinity. We are unable 
to take this continuum limit, for we must maintain the value of the excitation gap of the 
unpaired down electrons A^, ~ V^/D: this means that we must maintain a finite band 
electron cutoff. (In other words, the number of electrons per local moment ~ N{0)D must 
remain finite.) Since our mean field theory will not be accurate at frequency scales that 
are large compared with the Kondo temperature, we shall choose a frequency cutoff that 
is intermediate between the conduction electron bandwidth and the Kondo temperature 

Tk «A«D 
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With this choice, we are able to replace the density of states by its value at the Fermi 
surface: the energy integral is then carried out in the same fashion as section (4.), replacing 



(5.19) 



where the contour integral proceeds clockwise around the branch cut in the function Q[z] = 



In 



z-D 



The energy integrals can then be performed by closing the contour around the 



poles in the Green functions, which are located at 

^^ai'^) =UJ-^u + a Sgn(w)^A2 +//2 
Carrying out the contour integral in 2; = e then gives 



(« = ±) 



(5.20) 



aA 
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I ^Im{0[ea(c. + ^^)]}th[/3c^/2]- 



(5.21) 



(A2 +^2)2 

These two poles cancel one another's contributions, except in the gap region |c<;| < A^,, 
where the gapless excitation branch contributes a finite amount to the stiffness. Our final 
result for the London stiffness is then 



Q = 4e^p3 



4m 





0U1 


/ dwih 


Jo 


[ 2 J 



(A2 +yu2)2 



F2 



(5.22) 



where we have set ^ = By making the low temperature expansion 



th 



(5.23) 



the temperature dependence of the Meissner stiffness and penetration depth become 



1 



1 



2rp2 



/X, TT^T 



3A2 



Ai(T) (A«)2 

where 4:TrQ{T) = [Xi{T)]~'^ defines the London penetration depth and 

1 



F{x) = 2x' 



1 - 



Vl + 4a;2 



(5.24) 



32 



This variation of the penetration depth is similar to that expected for point-nodes in a 



conventional pairing scenario. In the special case of /i = 0, the Meissner stiffness is simply 



suppressed by a factor of Z = ^g/D compared with a conventional metal. Loosely speak- 
ing, we may consider this to be an effect of the condensation of heavy fermions, whose 
effective mass is enhanced by a factor m*/m ~ 1/Z, and whose rigidity is then depressed 
by the factor m/m* ~ Z. 

The staggered phase of the order parameter plays a critical role in developing this finite 
stiffness. This point is illustrated in Fig. 8. The "uniform" odd-frequency triplet state 
{Q = 0) is unstable and its energy may be monotonically reduced by twisting the order 
parameter until the stable minimum at Q — (tt, tt, tt) is reached. 

The "spin stiffness" pm = Pi, 2 for twisting the order parameter about the dj_ axes can 
be calculated in a similar fashion. When // = 0, the system is particle-hole symmetric, 
and Pm = P3 as given above. Like the superfluid stiffness ps, contributions to the stiffness 
come predominantly from the neutral excitation band inside the gap, though the formal 
expression for /j, ^ is more complicated, and shall not be given here. 

5. (c) Collective Modes 

To end our discussion on the long-wavelength properties, we should like to briefly 
mention the collective modes of the condensate. Let us generalize the effective action to 
incorporate the leading order time dependence of the pairing field 




(5.25) 



where we have set = V /D. The stiffness of the order parameter is thus finite, but 




(5.26) 
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where Xm denotes the magnetic, or spin susceptibihty, and x denotes the charge suscepti- 
bihty. An apphed vector potential is gauge equivalent to a rotation about the axis, and 
must be included with the kinetic terms to maintain the gauge invariance, as explained 
above. The spin and charge susceptibilities are given by x = Xm = P- In the absence of 
a coupling to the electromagnetic field, this action would give rise to a collective phase 
mode, and a spin-wave mode, with velocities 



^Ln ~ ise ~ 4 ( ^ ) (5-27) 



'Tk 

■'spin ^phase \ £) 

Of course, phase modes are gauged away as fluctuations of the electromagnetic field, con- 
verting the phase mode into a longitudinal plasmon mode as part of the Meissner effect, 
following the well-known Anderson-Higgs mechanism.^^ The spin-wave mode cannot be 
gauged away in this fashion and is unscreened, leading to gapless collective spin modes 
in the superconducting state that coexist with the superconductivity. From the velocity 
of the spin- wave excitations, we can see that these modes cross into the bottom of the 
quasiparticle continuum at a wave vector 



where a is the lattice parameter. This is much smaller than the size of the Brillouin zone so 
long-wavelength fluctuations of the order parameter will not lead to a dramatic reduction 
in its magnitude. On length scales shorter than 



and at frequencies greater than 



e~a^/^ (5.29) 



U! ~ Ag 



this system will behave much in the way of a single Kondo impurity. The development of 
coherence on longer length-scales provides a vital cut-off to the infrared fluctuations that 
destroy the condensate in a single impurity model. ^-"^'^^ 
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6. EflPect of coherence factors on local Response Functions 



The unusual nature of the coherence factors in the quasiparticle excitations have inter- 
esting consequences for the low frequency response of this system. Of particular interest 
here, are the local dynamical spin susceptibility and charge susceptibilities 

fOO 

Jo 



ILVt 



(6.1) 



These functions are directly related to the ultrasonic attenuation and the NMR relaxation 
rate, 1/Ti 

(6.2) 



Q!s(T) =Ai lim 



— {b) =X2 lim 
J 1 w— >0 



X+-('^) 



associated with the conduction electrons. Let us focus on contributions to these response 
functions derived from the gapless excitations in the "up" spin band of our toy model. 

The imaginary part of the local spin or charge response function of these excitations 
is given by 

X%u;) 4xL(a;) 



UJ 



ki, k2 



CO 



(6.3) 



ki' 



(The only component of the susceptibility matrix which couples to the low energy quasi 
particles is Xzz-) Since the coherence factors grow linearly in the energy 



The low energy form of this response function is given by 



(6.4) 



{...} 



TT j dEN%E+)Z{E+)N*(E-)Z(E-){...} 



f{E.) - f{E+) 



UJ 



{E+ + E-)^ 
F2 



(6.5) 
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where E± = E±oj/2. The density of states of the quasiparticles is N*{E) = {p/2)Z ^{E), 
thus at low temperatures and frequencies 

X=M/. = 4.LM/. = ^(^^^i^) (6.6) 

This quadratic temperature and frequency dependence of the local charge and spin re- 
sponses results from the unique energy dependence of the spin and charge matrix elements, 
and the neutrality and spinless character of excitations at the Fermi surface. In a more 
conventional superconductor, where gapless excitations carry charge and spin, these ma- 
trix elements would be unity at the Fermi energy, and this kind of quadratic behavior can 
only be produced by a linear density of states. A quadratic growth of the dynamic spin 
and charge susceptibility results in characteristic response of the NMR relaxation rate 
(Fig. 9). 

and a response of the ultrasonic attenuation rate. 

This property of the odd-frequency paired state is of particular interest, because power- 
law behavior of the above variety is observed in heavy fermion systems. Conventionally, 
it is ascribed to d-wave pairing and gaps vanishing along lines of the Fermi surface. Our 
odd-paired state offers the interesting alternative description of this as a matrix element 
effect. 

As a more specific illustration of this effect, let us examine the detailed temperature 
dependent NMR relaxation effect in our toy model. The complete local dynamical spin 
susceptibility has contributions from both the conduction electrons and the local spins 

Xii{oj) = -i [ dteM-i^t)e{t){[Si{t),Si{0)]); 

^ v-r- -1 (6.9) 

i 
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In general, the local moment spin operator mixes the gapless and gapful excitations and 
only the conduction electron spin operator in the z direction is able to create gapless exci- 
tations. Let us examine the corresponding conduction electron spin response in Matsubara 
frequency, given simply by 



Xlzi'^^n) = {Sz{'ioJn)Sz{-ioJn)) = ^ X] (pcaiii^n) Pcai-ii^n)) (6.10) 
We may now write 



where we have used the summations over momentum to replace the "up" conduction 
propagators by the corresponding local propagator (3.29) 

k 

where 

Ecu = sgn(a;)^A2 +//2. 

PT«(^) = (Vp)M^Ta(^ + ^<^)]- 
This spectral density can be used to carry out the Matsubara frequency convolution. 
We can consider the a = + term only, since it corresponds to the gapless excitations. The 
final expression for the imaginary part of the low-energy spin response is 

xi'zi^) _ V f ^^^^ j^x)- f{x-u)_f I 

(6.11) 



to oZi J CU 



{...} 



1 



ExEx—u) ' 

At low energies, the energy dependence of the magnetic relaxation is governed by the 
linear energy dependence of the spin coherence factors. The low energy imaginary part of 
the susceptibility is hence quadratic in the external frequency. 



u 48 \d) \~ 
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This in turn leads to an NMR relaxation rate that is proportional to 



1 



Trp^ / \2 /27rT\ 2 



TiT 48 \D/ V Tx 

All other parts of the spin response function lead to a gapped contribution to the 
magnetic relaxation rate. Furthermore, in the undoped case, the mid-gap coherence factors 
are identically zero, and the imaginary part of the spin and charge response functions are 
zero throughout the gap (Fig. 10). 

The key difference between the d-wave and odd frequency description, is that the 
former relies on a node in momentum space, whereas the latter relies on a node in frequency 
space. For this reason, we expect that odd-frequency pairing is rather insensitive to elastic 
scattering. This point can be illustrated in the following general way. Let us consider a 
general odd-frequency paired state, where the pairing self-energy takes the form 

/ 

X, X ) ^Q{x)—V{x)5xx' 

^ (6.12) 

V =^ [3(1) - dab{x)a^ ® T^] ; 

here the amplitude VqC,{x) of the resonant scattering and the orientation (^^6 = 
{x)aaTi,Z{x) of the triplet order parameter may be site dependent. The component 
of the conduction electron states which does not directly couple to the resonant pairing, 
'^oi.x) may be projected out of the conduction electron spinor as follows 

^o{x) =p{x)^{x) 



p{x) =1 - E(x) = ^ (1) + dab{x)a'' ® 



(6.13) 



This component experiences an indirect effect of the resonant pairing through mixing with 
the directly scattered components. 

Consider a general conduction electron band with disorder, described by the Hamilto- 
nian 
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Let us project out the parts of the conduction electron Hamiltonian that couple directly, 
or indirectly to the resonant scattering, writing 



h 



XX' 



t , 



a 



x.x' 



p{x) 



V{x) 



H{x, x) 



p{x) 



V{x) 



(6.15) 



The conduction electron Green function can then be written 

-1 



uj-n- -^cv 

UJ - 



^^xx' = C(a;)^a;,x') 



(6.16) 



where all subscripts have been omitted. The projected component of the conduction elec- 
tron propagator for those states that do not directly couple to the resonant scattering is 
then given by 

G-l(a;) =a;-/i-a(a;-ii'- ^C) (6-17) 

Though these states do not directly couple to the resonant scattering, they couple indirectly 
because of the off-diagonal terms in the Hamiltonian that mix them with the resonantly 
scattered states. At low energies, the resonant scattering dominates the "self-energy" 
correction in the propagator, which then becomes 

G^l(a;) =Zr^u-h 



1 + 



(6.18) 



Notice that the effect of the resonant scattering is to introduce a wavefunction renormal- 
ization into the propagator. Low energy eigenstates are set by the determinantal equation 



^et[G-l] = 



Clearly then, if there are zero energy eigenstates of the projected Hamiltonian 



(6.19) 



h{x,x^)iyix') = ^ 



(6.20) 



then these will give rise to zeroes in this determinant. In other words: the projective 
character of the resonant scattering means that the indirectly coupled zero energy states 
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form zero energy excitations of the complete Hamiltonian. Suppose we define the Majorana 
conduction electron states 

ax'f = J2M^)U^) (6.21) 

X 

Then their propagator will be given by 



(a(a;)at(a;)) = xG^i^)^X = — (6-22) 
where the pole strength Zx is 

l + ^x^-p^ix (6.23) 



Thus the zero energy eigenstates of the complete Hamiltonian will have the form 



= ^J2MxKx(x) + ■ ■ ■ (6.24) 



where the residual part of strength ^/(l — Zx) is carried entirely by the Majorana spin 
fermions. In general then, the off-diagonal coupling between the indirectly, and directly 
scattered states leads to a reduction in the conduction character of the Majorana zero 
modes, and a corresponding enhancement of the density of gapless excitations. 

Ar*(0) = (0 {Z-')x (6.25) 

The important point however, is that despite these effects, the gapless excitations remain 
Majorana fermions: the spin and charge operators are completely off-diagonal at the Fermi 
surface, and coherence factors must consistently vanish in this region. 

To provide a specific example, consider the generalization of our toy model with a 
random chemical potential 

He *"^k%^k ^Hx)^i{x)T^^{x) (6.26) 



k ^ 



2 ^ k k^k 2 
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The chemical potential term can be identified as the off'-diagonal coupling, whereas the 
kinetic energy term commutes with the projection. The Z-factor for a gapless plane wave 

Z.^=(l + {lj?{x))/V^) (6.27) 
which gives rise to an enhancement of the gapless density of states given by 

N\Q) = P-(l + {^\x))/V^) (6.28) 

where for weak scattering, we have ignored the fluctuations in the strength of the resonant 
scattering potential that will be induced by the disorder. Thus we see that disorder eu- 
hances the density of states but sustains the electric and magnetic neutrality of the Fermi 
surface. 

7. Interplay with magnetism 

In this section we discuss some of the magnetic aspects of the odd-frequency state. 
Even in our toy model, where we have not included any detailed effects of band-structure 
or anisotropy, there are a variety of locally stable phases where the order parameter is 
commensurately staggered. Quite generally, as we now show, the odd-frequency state will 
develop an ordered magnetic moment, aligned parallel to the vector, giving rise to 
coexistence of antiferromagnetism and superconductivity. 

7(a) Pair- Spin Correlations 

We begin our discussion by returning to the simplest example, where the (fi vector is 
uniformly oriented, giving rise to a state with ferromagnetic correlations. We will general- 
ize our discussions to a more realistic antiferromagnetically ordered case at the end of the 
section. When the conduction band is half filled, and hence completely particle-hole sym- 
metric, the odd-frequency paired state is magnetically isotropic and both static magnetic 
order and static pair correlations are absent (S) = 0, {tj) = at half-filling. In this state, 
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there is long range order with an order parameter 

{S''{x)r\x)) = 2^ [d,{x))\ 

where S{x) = Sc{x) + Sj:{x) is the total moment at site x, and t{x) is the conduction 
electron isospin at site x. Odd frequency triplet pairing is thus seen to strongly couple 
spin and pair correlations. Clearly however, these cross-correlations induce anomalous 
response functions, coupling the development of charge correlations to the application of a 
magnetic field or the development of magnetic correlations to the application of a chemical 
potential or pairing field. To study this effect, we introduce the spin-charge susceptibility 

where the static susceptibility x(0) is of particular interest. In the vicinity of half filling, 
the state is completely isotropic in spin and isospin space, and we then expect 

X"^ = Xod^^ (6.29) 

Thus, once the system is doped (// 7^ 0), the presence of a cross-correlation between the 
charge and spin degrees of freedom leads to the development of a magnetic moment and 
strong anisotropy 

We may calculate Xo very simply as follows. Let us choose cfi = z. The coupling of the 
magnetic field is given by 

Hb^-J2^- - \vj X flj] (6.30) 

j 

Thus the total magnetic moment in the z direction per site is 

1 (6.31) 

= (2(M - M) + (V2 - rirji)) 
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here rica, ( cr =1, i) is the number of "up" or "down" conduction electrons per site and 
following the notation of (3.17), n,j| = jy^^ ^k^kj^^kj, number of "down" Majo- 

rana fermions per site. Now note that the down Majorana electrons are hybridized with 
the down conduction electrons to produce a completely filled hybridized band, with one 
electron per unit cell. Thus in the ground-state, {rirj + n^^) = 1 per site, so 

Mz = linc^ + n,^) - ^ (6.32) 

per site. Since the local moments leave the conduction electron density essentially un- 
changed, (nc| + rici - 1) = 2p/i, giving 

Mz = pii (6.33) 

and 

Xo^P + 0{Tk/D) (6.34) 

A departure from particle hole symmetry Sfi thus generates an ordered (ferromagnetic) 
moment of strength 

M(x) = 5iipd^(x) (6.35) 

In a similar fashion, an application of a magnetic field will influence the charge and 
pair correlations: a field along the direction will develop a charge density, 

p{x) = pBz{x) 

More remarkably, a transverse magnetic field B — Bj_d}, (/ = 1,2) will induce a conven- 
tional pairing field of magnitude 



This leads to a possible field dependent Josephson coupling with conventional supercon- 
ductors. 
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Even in this ferromagnetically ordered state, there are strong antiferromagnetic corre- 
lations. To illustrate this point, consider the case where n — 0. Here, the main contribution 
to the magnetic susceptibility is provided by the local moments, and may be calculated 
from the spatial correlations of the Majorana fields. Low frequency properties of the local 
moments are determined by excitations across the indirect gap at wave vectors q — ±Q/2. 
In the vicinity of the gap, the Majorana fermions can be expanded as shown in (3.5) and 
(3.6). In momentum space, the Majorana propagators have the form 

D{k) 



- c^) (6.36) 



[ioJn{ioJn - e^) - V^] 
so that 

ImD(k, u! — i6) = ^2 — 2'^S{u! — a;^) (6.37) 
The spin correlations are determined from the product of two Majorana propagators: 

{S^iq)S\-q)) = S^\iq) = S^^^ J] Diq/2 + k)D{k - q/2) (6.38) 

At low frequencies, this is dominated by governed by excitations in the vicinity of the 
indirect gap, the quasiparticle spectrum is parabolic ujq = ±A^ + ; where m* = 



D 

J-K 



m, and m is the electron mass at the band-edge. Using this parabolic approximation 



to carry out the momentum space integrals gives 

y2 



Imx{q,u; + iS) 



y2 + (a;/2)2 



(A, =2A, + i^^) 
^ q 4m* ^ 



Ju - A^eiu - Ag) 

(6.39) 



In other words, independently of the ordered magnetic moments, there is a large amount 
of spin-fluctuation spectral weight above the superconducting gap at the antiferromag- 
netic zone vector. This general feature survives when we come to consider more general, 
staggered configurations of d^{x). 
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Let us now consider the possibility of more general, staggered configurations of cfi. 
Take the more general ansatz for the mean field order parameter 



di{x) \ / Jicos[Q-:^ \ (6.40) 
M{x) = I cUx) ] ■ = [ d2 cos[(g + P)-x] . 

— * — ♦ 

where P and Q are commensurate vectors. 

As before, we can redefine the conduction electron states to take account of the stag- 
gered order. 



The conduction electron Hamiltonian can then be written 



(6.41) 



Hc= Yl *^k[(e[k-Q/2r3-P/2^i]-/.)T3]*j^ (6.42) 



keiB.z. 



where the kinetic energy term can be expanded as 
[e(k-g/2T3-P/2ai)-/x]T3 



--\ E {(T3 + «)(l + /?^i)6[k-ag/2-/?Q/2]}-/.T3 



(6.43) 



From the discussion of section 3, we know that gapless excitations will develop on the 
"Fermi surfaces" described by 



J J2 ae[k -aQ/2- (3P/2] = (6.44) 



i 4 



To make our example more specific, consider the case where 

e^ = - 2t[cx + cy + Cz] {ci = cos[ki],l = 1, 2, 3) 

Q=(7r,7r,7r), P = (tt, 0, 0) 
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(6.45) 



corresponding to a staggered cfi vector in the x direction. In this case, 




2t[Sy + Sz] 



{si = sm[ki],l = l,2,S) 



ei =2tc. 



k 



'X 



(6.46) 





The gapless modes he on a tube with a square cross-section [sy + s^] = and the spectrum 



is given by 



Det G-^{uj)G-^{uj) - [el]^! = 
G-l=[a;-ee-2A^ + /XT3] 



(6.47) 



After a short calculation the corresponding mass renormalization of the gapless quasipar- 
ticles is found to be 



clearly showing the development of local minima at the commensurate points in the Bril- 
louin zone. At each of these points the odd-frequency state will develop a staggered mag- 

— * 

netization with wave vector P and approximate magnitude M ~ pfi. In practice, we would 
argue that the small differences in energy between these different commensurate states will 
depend on several factors that are not included in the toy model. It is interesting to note 

— * 

that in two dimensions, the toy model predicts that a P = (tt, 0) state is more stable than 
the P = (tt, tt) state. The important point however, is that the mean field energy of these 
magnetic phases is close in energy to the uniform state and furthermore, is locally sta- 
ble. We may conclude that this type of odd frequency pairing can homogeneously coexist 
with antiferromagnetism. The local moments participate in both the spin and the pair 
condensate. 




(6.48) 



In Fig. 11, we show the mean-field ground-state energy as a function of P = {Px, Py, Pz), 
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8. Critique and Discussion 

Our paper has presented a "toy" realization of odd frequency pairing, with the aim 
of elucidating its key properties. In this section, we discuss the odd frequency state in a 
more general setting and examine the possibility that this kind of paired state might be 
applicable to heavy fermions. 

One of the most dramatic features of the theory is the projective character of the 
resonant pairing self-energy: 

S(a;) =—V 

^ (6.49) 

Can we understand this feature in a more general context, outside the restrictive realm of 
the Kondo model, and our Majorana treatment? Let us consider the possible extension 
of our odd paired state within an Anderson model for heavy fermions, with an on-site 
repulsion term 

Hi = ]^U{n-lf 

at each magnetic site. Highly correlated states minimize this on-site interaction energy, 
tending to produce local moment states where n = 1. It is quite useful to examine this 
constraint in terms of the correlations between the charge (isospin) and spin of the localized 
states. If we expand the localized f-electron in terms of its four real components 

where, z = ^ ^| ^ ^ spinor, then the interaction may be written as a symmetric 
product of all four fields. '^'^ 

\ (6.51) 
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On a lattice, the interaction between the electrons can then be written 

= -Uj2{flfkfifl}^^o+n^+n,+ns (6-52) 

K 

where we use the Fourier transformed operators 

r(^) = -^^ J dre-'^-^fiX), K-X = k.Rj-uJnr. (6.53) 

and it is understood that all time coordinates are ultimately to be time-ordered. 

A novel way to reduce the on-site correlation energy is to develop a correlated state 
where certain Majorana components of the f-state are absent. Remarkably, the operator 
that projects out the zeroth component at wave vector a 

f{K) = ^[z'ff{K) + fH-K)z\ (6.54) 

is a one particle operator 
where d^b = \z^aa ® t^Z^ 

z = ( . 



z 

and 



z 



are the Balian Werthammer four spinors for the f-state and z spinor. The residual "vector" 
(1, 2, 3) components of the f-states are projected out by the one particle operator 

i=i,3 

The Majorana character of this operator implies that it is asymmetric 

A = -V-n (6.55) 
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The self-energy term that will selectively decouple Majorana components of the f- 
electrons will have the general form 



(6.56) 



Since Vk is an odd function of k, it follows that 

A(k) = -A(-k). (6.57) 

If the physics is local in time, then the frequency dependence of A{k,) can be dropped, 
leading to p-wave triplet pairing. However, if the momentum dependence of A(k) is even, 
the frequency dependence is automatically odd, leading to odd frequency pairing. In the 
simplest s-wave version of this pairing, the physics is spatially local, so that A(k) = A(a;). 
This establishes an intimate connection between the projection of Majorana degrees of 
freedom from the ground-state and the development of nodes in the wavefunction: when 
the physics is local, this projection results in node in time, and the development of odd- 
frequency pairing. 

A general spectral decomposition of A(a;) will always contain a zero frequency pole 

a; y 7ra;2-z.2 ^^^^^^ 

-Im[A{uj - i6)] ^Z6{u;) + A{u;) 

TT 

This pole is a unique feature of odd-frequency pairing: it suppresses the "vector" com- 
ponents of the f-electron from the low energy excitations, gaining correlation energy and 
decoupling a band of gapless singlet excitations. Our simple mean field theory can be 
viewed as a dominant pole approximation to the pairing field. A pressing need for the 
near future is to show that such general constructions can lead to stable Eliashberg-type 
treatments of more general models, such as the finite U Anderson lattice. 
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In its current form, our prototype for odd frequency pairing is too simplistic to account 
for details of heavy fermion behavior. We should like to list some important issues which 
need to be addressed in future developments: 

• Magnetism. The toy model has shown that odd-frequency pairing has a propensity 
to coexist with magnetism. A more realistic model will need to take explicit account 
of the RKKY interactions and their role in establishing the detailed superconducting 
order. 

• Normal Phase. The normal phase of heavy fermion superconductors, with its pro- 
fusion of Fermi liquid features, does not appear in the toy model. The Majorana 
formalism may be a poor starting point to recover the normal phase properties, and 
this suggests that we should seek a way to obtain the odd-frequency paired state within 
perturbation theory for the finite U Anderson model, or perhaps the large N approach 
to the heavy fermion problem. 

• Anisotropy. Measurements of the gap and ultrasound absorption^^'^^ in U Pt^ show 
the presence of anisotropy in the order parameter and have been traditionally inter- 
preted within a d-wave pairing scenario. These results do not reveal the temporal parity 
of the paired state, but tend to reinforce the conclusion that momentum anisotropy 
can not be ignored in a more advanced version of the model. Indeed, there is no reason 
not to contemplate the possibility of odd frequency d-wave pairing. 

• Power Laws. Power laws in the temperature dependence of the specific heat and 
NMR relaxation rate of heavy Fermion compounds develop much closer to Tc than 
any simple mean field theory can account for. One possibility is that dynamic pair- 
breaking effects have suppressed Tc significantly below the gap. Odd frequency pairing 
accounts for the finite linear specific heats in heavy fermion superconductors in terms 
of a band of excitations with vanishing coherence factors. At present, the toy model is 
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unable to account for the term in Cy that is also seen. Fluctuation effects need to 
be examined carefully. 

We should like to spend a moment discussing the long-wavelength properties of odd- 
frequency triplet paired states. The intimate relation between spin and pair degrees of 
freedom in this kind of state leads to rather interesting consequences in the Landau- 
Ginzburg theory. Suppose one considers the simple long-wavelength action discussed in 
section (5) 



(fix 



pA^nf + Ps{i^3 - j^Af + r_J _ - y d^xid^ ■ zY (6.59) 

then the anisotropy plays a vital role in establishing the topological stability of persistent 
currents. Unlike a conventional superconductor, the supercurrent is linked to the spin order 
and involves all three Euler angles of the order parameter. To see this, it is instructive to 
consider a loop of superconductor of length L, threaded by a solenoid. The supercurrent 
around the loop is given by 

2e _ 2e 

^ ^ (6.60) 

0)3 =V0 + cos OVip 

where (0, 9, 0) are the Euler angles defining the orientation of the triad d^^^. Notice that 
unlike a conventional superconductor, the supercurrent involves both the U(l) phase (j) 
and the orientation of the magnetic vector d, defined by (OjI/j). For a conventional super- 
conductor the total phase change around the loop is a topological invariant 



2e -9 



A0 



j dr-V^ = 27rn (6.61) 



that is unchanged upon application of a fiux through the solenoid, leading to a linear 
relation js = ~^^Ps-§- between the enclosed flux $ and the supercurrent density js- In 
this superconductor, the analogous integral around the current loop is 

A(f) = J dr- V(f) + cos eWip = 27rn + n (6.62) 
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The second term 

n = J dr- cos ^VV'] = J dS- {eab^an x V^n) (6.63) 

is the sohd angle subtended by the h vector around the loop: this is not an invariant, and 
can change by multiples of 47r to relax the current. Unlike the U{1) superconductor, the 
only stable vortex configuration involves a net phase change of 27r around the loop: this is 
the so called "Z2" vortex of an S0(3) order parameter, and it has the property that two 
such vortices can be adiabatically deformed back to the vacuum (Fig. 12). Each Z2 vortex 
pair reduces the effective flux through the solenoid by 2 flux quanta, thus the non-linear 
current will be given by 



2e 

where the number of Z2 vortex pairs is 



27r— -47rn^J (6.64) 



nz = Int 



$ 1 
+ - 



(6.65) 



-2^0 2. 

where Int(a;) denotes the largest integer smaller than x. In this way, the current density 
around the loop can never exceed 

• _ (27r)2 

and will never become macroscopic. Thus, without anisotropy, the critical current for the 
odd-frequency state is zero. Furthermore, if the fliix through the solenoid is $ = J V{t)dt., 
then the response to an oscillatory electric field will not occur at the driving frequency. This 
will eliminate the low frequency linear Meissner response to macroscopic fields, removing 
the pole in the optical conductivity and producing an apparent violation of the linear 
response optical sum rule. This type of behavior is most likely to occur in the vicinity of 
particle-hole symmetry, and suggests that this half filled state will more closely resemble an 
insulator, rather than a superconductor. This may be an interesting way of thinking about 
Kondo insulators^, where an anomalous reduction in the low frequency oscillator strength 
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of the optical conductivity has recently been reported. " As anisotropy is increased, a 
macroscopic Free energy barrier will have to be crossed in order to add pairs of Z2 vortices 
to the superconducting state, restoring the linear Meissner response. This will lead to a 
non-trivial dependence of the critical current on anisotropy. 

Pending further theoretical work, it appears that there may be some useful exper- 
iments on heavy fermion superconductors that could help to compare the d-wave and 
odd-frequency scenarios. A key issue is to verify the relation between the gapless excita- 
tions and the NMR relaxation rate, most notably to confirm the presence or absence of a 
Korringa term in the relaxation rate in severely gapless heavy fermion superconductors. 
Another area of fruitful investigation concerns the field dependence of the proximity efi^ect. 
Negative proximity effects have been observed between U Bei^ and To superconductors.^^ 
If the symmetry of heavy fermion superconductors has a different temporal parity, then we 
expect the application of a magnetic field to enhance the coupling between the two order 
parameters, leading to a strong reduction of the negative proximity effect in a field, and a 
strong field dependence of the Josephson current. 

In conclusion, we have presented a stable realization of odd frequency triplet pairing 
in a Kondo lattice model for heavy fermions. Under rather general conditions, the odd 
frequency state that forms has a gapless singlet mode of quasipar tides. Spin and charge 
coherence factors for these quasiparticles grow linearly in their energy. Our pairing hy- 
pothesis provides an alternative explanation of various power laws in heavy fermions in 
terms of a vanishing of coherence factors at the Fermi energy, rather than a vanishing of 
the density of states. We have conjectured that this may explain the absence of a Korringa 
law in the NMR, even when the superconductor is highly gapless. Odd- frequency triplet 
superconductivity appears to be able to coexist with magnetism, and in our simple toy 
model, ferromagnetic order coexists with the pairing. We think our results are encour- 
aging enough to prompt efforts to develop a description of odd-frequency pairing within 
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more general models, and to consider seriously the possibility that this provides a viable 
alternative to the d-wave pairing hypothesis in heavy fermion superconductors. 
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APPENDIX A. Majorana Representation of Spins 

In this section, we present a derivation of the Majorana representation that provides a 
link with the Abrikosov fermion representation and illustrates how the constraint is avoided 
by the uniform replication of the spin Hilbert space. We begin by noting that for any two 
component electron spinor there are two operators of interest: the "spin" 



and the "isospin" 



ZAap 



2 



a/3 



where we have introduced the Nambu spinor 



/ = I /t\ 1 • (A.3) 
These operators are independent [sa,T^ = and each satisfy an SU{2) algebra 

[si, Sj] = ie'^'^Sk, [Ti, Tj] = ie'^'^Tk, (A.4) 

In the subspace where the spin is finite, the isospin is zero, and vice versa. The sum of 
both operators 

S = s + T, (A.5) 
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satisfies an SU{2) algebra, and is either equal to the "spin" or "isospin", depending on 
which component of the Fock space is projected. For any interacting system of electrons 
containing N local moments, we may write the partition function as a constrained trace 



(qj = s, t) 



with 



4- = 



a(3 



t. 



ja 



a(3 



(A.6) 



(A.7) 



where the projection operator projects the "spin" or "isospin" component of the Fock 



space at site j 



PJ ={nj - ly 



(A.8) 



n 



There are then 2^ ways of choosing the projection operators: each choice projects a 
replica of the spin Hilbert space with precisely the same partition function. Summing over 
all replicas we can write 



(A.9) 



Qj=S,T 



The sum over all 2^ projectors is the identity operator 



(A.IO) 



and hence the replicated partition function can be written as an unconstrained trace, with 
each local moment represented as a sum of the Pauli spin and isospin. 



Z = 4fTr 
2N 



(A.ll) 



We now demonstrate that the combined operator S = s + r depends only on three Ma- 
jorana components of the f-electron. Suppose we decompose the complex Fermi operators 
into their real and imaginary Majorana components as follows 



(A.12) 
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In terms of these components, the "spin" and "isospin" operators are 



1 



m 



(A.13) 



The sum of these two is then 



+ T • 



_]iJlm I m 
2 'j 'j 



(A.14) 



which is precisely our Majorana representation. With this choice of z, the zeroth com- 
ponent of the Majorana fermions at each site does not enter into the Hamiltonian. This 
component can therefore be exphcitly traced out of the partition function. Formally, this 
may be done by pairing the zeroth Majorana fermions throughout the lattice 



(A = l,2...iV/2) 



(A.15) 



where each site / belongs to one pair: / e {(^AiJa)'-^ = li-^}- The set of N/2 complex 
fermions are independent and form a completely decoupled zero energy Fock space of 
dimension Hence 



Z = 



1 



2N/2 



Tr 



(A.16) 



where the remaining unconstrained trace is over the I = 1, 2, 3 components of each Majo- 
rana fermion, and the other real electron states of the system. 

The overcompleteness of our representation is closely related to a residual discrete local 
Z2 symmetry of the Majorana spin representation under the transformation 



(A.17) 



In this respect, the Majorana representation is similar to the pseudo fermion representa- 
tion. However, in this case, the canonical and grand-canonical ensembles have precisely 
the same partition function, up to a simple normalization factor. In the pseudofermion 
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representation, the Gibbs partition function for each conserved subspace is not constant 
(for if nj = or nj = 2 there is no spin at site j), and the projection of the unwanted 
spaces is an unavoidable necessity. 

Majorana fermions are easily treated in a momentum space representation. The Fourier 
transformed operators 

i 

Since the original Majorana fermions are real r]j = rf^\ their complex Fourier transforms 

satisfy ryS = ■n'^t., forming a set of independent complex fermions that span half of the 
k — k 

Brillouin Zone 

rf^ I = 5^^^ 3^ ^ k, k e half the Brillouin Zone. (A. 19) 

The inverse transformation can be written 

kel/2B.Z. 

The corresponding Lagrangian for the Majoranas is then 

^= E %^-drX + 'n (A.21) 
keiB. z. 

or in terms of the original site representation 

C=\Y.^i-drm + U. (A.22) 

i 

Note that for each momentum fc, we can choose either ry^ or r]_^ — ryt- as the independent 
destruction operator. This has an important consequence for broken symmetry solutions, 
for there are equivalent ways of making the choice of the vacuum state: by making one 
particular choice, the normalization constant in front of the partition function is absorbed. 

APPENDIX B. Some simple examples 
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In this Appendix, we illustrate the use of the Majorana fermion representation of spin 
1/2, by means of some specific examples. 

Consider first the Heisenberg Model for 2 spins 1/2 (where we take the exchange 
coupling to be 1), written in terms of the Majorana fermions 

H = Si-S2 = ~[{m-m? + l]. (B.i) 

One can now define 3 complex fermions, by taking appropriate linear combinations of the 
6 Majorana fermions 

V (B.2) 

These operators satisfy the usual fermionic anticommutation algebra 

{f\f^ = Sij {iJ = l,S) (B.3) 

and act on a Hilbert space of dimension 2"^ = 8. As pointed out in the first section, the 
dimensionality of the original space has been increased by a factor of 2^-^/'^^ = 2, where 
jV = 2 is the number of spins. 

By using 

-im ■m = P ■ / - ^, (B.4) 
one can now write the Hamiltonian in terms of the f-operators 

8 2^"^ •' 2' 

The spectrum can now be easily worked out 

Eo = -3/4, 

(B.5) 

El = 1/4. 

The first level is doubly degenerate and the second one is 6-fold degenerate. The exact 
eigenenergies are correctly obtained, as expected. Besides, the singlet ground state and 
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the triplet excited state are replicated by the same factor of 2. This additional degener- 
acy can be traced back to the invariance of the Majorana representation with respect to 
the Zi transformations r/^ —ffi, which is reflected in a particle- hole symmetry of the 
Hamiltonian f ^ ■ 

A mean-field treatment of this model can be performed by the following decoupling 
procedure 

Hmf = 1 + iVim ■m) + ^ = l- viff . / - ^) + ^. (B.6) 

A static order parameter V breaks the aforementioned Z2 symmetry and there are two 
stable solutions related to each other by the transformation 

V^-V, f^ff. (B.7) 

If V is positive, the ground state corresponds to / ' • /— ^ — V — ^ and its energy is —3/4, 
which is the exact value. 

Consider now the one-dimensional XY model, which can be solved exactly through a 
mapping to a free fermion model. This mapping is conventionally performed by means of 
a Jordan- Wigner transformation. It will now be shown that an analogous mapping can be 
achieved through the Majorana representation of spins. 

The ID XY Model Hamiltonian is given by 

HxY = E i^'^^x-^' + ^i^i"') = 2 E {^+^-' + ^-^+') ■ 

i=l i=l 

In terms of the Majorana fermions, one can write 



SX^Sl, + iSl = rjM+iv'2), 



(B.9) 



si^si- isi = ivl - ivh)vl- 

At each site one can define a complex fermion 
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(B.IO) 



satisfying 

{c\ c^t} = 5i,; {c\ 4) = 0; {c'\ 4) = 0. (B.ll) 

The Hamiltonian can now be written as 

M 

HxY = E (^'"^44+'^'+' + c%vt^^+^^) . (B.12) 
i=l 

Let W be an operator acting on site i 

U' = Pi + V24pi (B.13) 

where 

= 1 - c'^c\ Pi = c'^c\ (B.14) 

are projectors onto the vacant and occupied states of site i. The operator C/* is both 
hermitian and unitary 

C/^t = P^ + V2Plr]l = U\ 

(B.15) 

Ijilji] = ui]ui = + 2r/f = + Pf = 1. 
One can now easily prove that 

^2 (B.16) 

The canonical transformation generated by transforms away the third Majorana com- 
ponent. 

Let U be the ordered product 

N 

U=WU' (B.17) 

i=l 

Under the action of [/, the transformed operator now acquires a non-local phase factor 

3=1 k=i j=l \^^/ k=i-l 



(B.18) 
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where and, in the last step, use has been made of 

t/i(c^t)c/i = ^\{pi _ V24p(){pi + V24p() 

(B.19) 

= c^t(_i)% ^i^jy 
Using relation (B.18) and its complex conjugate, one can transform the Hamiltonian 

into 



i=l 

1 ^ (B 20) 

i=l 

i=l 

which is the usual free fermion expression obtained by the Jordan- Wigner transformation. 
The third Majorana components have been transformed out of the problem. Tracing over 
these variables will cancel out the overall factor of 2-^/^ and one is left with a free fermion 
theory. 
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FIGURE CAPTIONS 



Fig. 1. Illustrating the elementary "Z2" vortex for a charge e spinor order parameter. 
Around the vortex, the phase change of the order parameter is tt. The supercurrent around 
this vortex can be removed by introducing a magnetic fliix for which ^ = ^ J A.dl = tt, 
so the flux quantum for a charge e spinor is $ = ^ = ^1, which corresponds to the flux 
quantum of a charge 2e complex scalar order parameter. 

Fig. 2. (i) Bare Majorana and conduction propagators, (ii) Interaction between local 
moment and conduction electrons. 

Fig. 3. Diagrammatic illustration of the pairing equations showing (i) the spinor vertex 
between conduction and Majorana spins (ii) the self consistent equation for the conven- 
tional and anomalous conduction electron propagators. 

Fig. 4. Schematic illustration of the quasiparticle excitation spectrum. The gapful spin 
excitations are separated from the gapless band of neutral singlet Majorana excitations by 
half the Brillouin zone. 

Fig. 5. Quasiparticle spectrum within mean field theory for /i/D = 1/6. Bold line: 
gapped "up" excitations; dashed line gapped "down" excitations; dotted line: neutral 
singlet Majorana band. Inset, density of states for up electron bands. 

Fig. 6. Illustrating the distinction between (i) the Hartree/RPA decoupling scheme used 
and (ii) a Hartree-Fock decoupling procedure. The "Fock" part of the vertex equation 
is absent from the Hartree approximation, but is reincorporated as a leading term in the 
RPA. 

Fig. 7. Temperature dependence of the specific heat, calculated at low temperatures for 
a variety of ///-D. Note that as n/D increases, the gapless band becomes "heavy". 

Fig. 8. Mean field free energy of the odd-frequency paired state plotted as a function of 
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the "twist" wave vector of the triplet pairing field dc{R) = e'^'^'^ddO). The uniform state 
{Q = 0) is unstable with respect to the state with a staggered phase {Q = n). 

Fig. 9. Conduction electron NMR relaxation rate, computed using our toy model for 
the same sequence of of /i/D values shown in Fig. 7. Even though the model predicts 
a linear specific heat (Fig. 7.), the coherence factors give rise to an NMR relaxation rate 
normally be associated with lines, rather than surfaces of gapless excitations. There is no 
Hebel-Slichter peak below the mean-field transition temperature. 

Fig. 10. Local dynamic spin susceptibility of the odd frequency state for various values 
of ji/D. Inset details mid-gap response that grows quadratically with the frequency, due 
to linear spin coherence factors. 

— * 

Fig. 11. Dependence of mean field free energy on the magnetic wave vector P in three 
dimensions. In three dimensions, the state at P = (tt, 0, 0) is locally stable. 

Fig. 12. Illustrating the energy dependence on the flux through a ring of odd frequency 
triplet superconductor (no spin anisotropy). The phase displayed here refers to the 
phase of the charge 2e composite order parameter gM.£{x) = {Ta(x)S^(x)). When the 
flux through the loop exceeds one flux quantum, the system can relax the energy and 
supercurrent by rotating the axes of the order parameter into the third dimension, creating 
two Z2 antivortices. This reduces the effective flux through the ring by 2 flux quanta, 
changing the sign of the supercurrent J = —dE/dcj) and producing a saw-tooth dependence 
of current on flux. The ground-state energy is a periodic, rather than a quadratic function 
of the applied flux. A macroscopic Meissner current can not develop in response to the 
flux threading the loop unless spin anisotropy is added to prevent the rotation of the order 
parameter into the third dimension. 
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